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Abstract We present the latest results from a multi- 
epoch timing and spectral study of the Transient Anoma- 
lous X-ray Pulsar XTE J1810-197. We have acquired 
seven observations of this pulsar with the Newton X-ray 
Multi-mirror Mission (XMM-Newton) over the course of 
two and a half years, to follow the spectral evolution 
as the source fades from outburst. The spectrum is ar- 
guably best characterized by a two-temperature black- 
body whose luminosities are decreasing exponentially with 
Ti = 870 d and T2 = 280 d, respectively. The tempera- 
tures of these components are currently cooling at a rate 
of 22% per year from a nearly constant value recorded 
at earlier epochs of kTi = 0.25 keV and kT2 = 0.67 keV, 
respectively. The new data show that the temperature Ti 
and luminosity of that component have nearly returned 
to their historic quiescent levels and that its pulsed frac- 
tion, which has steadily decreased with time, is now con- 
sistent with the previous lack of detected pulsations in 
quiescence. We also summarize the detections of radio 
emission from XTE J1810~197, the first confirmed for 
any AXP. We consider possible models for the emission 
geometry and mechanisms of XTE J1810-197. 
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1 Introduction 

Neutron star (NS) astronomy has been recently invigo- 
rated by the identification of a new class of magnetically 
dominated emitters. Known as anomalous X-ray pulsars 
(AXPs) and soft gamma-ray repeaters (SGRs), these ob- 
jects are apparently young, isolated neutron stars (NSs), 
whose properties differ markedly from those of the Crab 
pulsar, previously co nsidered prototypical o f the young 
NSs (for a review see 'Mere ehetti et al.ll200l . All AXPs 
and at least 3 of the 4 SGRs are identified as relatively 
slow (5 — 12 s) pulsars. Many are located at the cen- 
ters of recognized, young supernova remnants (SNRs), 
directly associating them with their supernova explo- 
sions. These objects emit predominantly at X-ray en- 
ergies and are distinguished by their characteristic spec- 
tral signature. This radiation cannot be accounted for by 
rotational energy losses alone, as for the radio pulsars, 
but is most likely powered by the decay of an enormous 
magnetic field characterized by a dipole with ^ 4.4 x 
10^^ G, at the pole. Collectively, these isolated NSs are 
understood within the context of the magnetar theory 
l)Duncan fc Thompsonll992HThompson fc Puncanl 1996(1 . 

A unique pulsar has been discovered whose remark- 
able properties offer great promise for deciphering the 
emission mechanism(s) of magnetars. XTE J1810— 197 
is a 5.54 s X-ray pulsar whose measured and derived 
physical parameters are fully characteristic of an AXP, 
but expresses b ehavior not previous ly associated with 
any such object l|Gotthelf et alJl2004l hereafter Paper I). 
XTE J1810-197 is a transient AXP (TAX P) - it was dis- 
cover ed during a bright impulsive outburst l|lbrahim et alJ 
I2004D that is still fading steadily. Even more surprising is 
the discovery of highly variable radio emission, providing 
the first confirmed example of radio flux from an AXP 
ySalpern et al. 2005b), and the sub sequent detection o f 
radio pulsations at the X-ray period l|Camilo et al ]200l). 

The outburst that resulted in the detection of the 
transient AXP occurre d sometime between 2002 Novem- 
ber and 2003 January l|lbrahim et alJl2'004|) . Since then. 
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over the course of a year, regular scans of the region 
with RXTE recorded an exponential flux decay with 
a time-constant of 269 ± 25 days from a maximum of 
F(2 — 10 keV) « 8 X 10"^-^ erg cm~^ s^^. In comparison, 
the previous average quiescent flux, with its softer spec- 
trum, gives F(0.5 — 10 keV) « 5.5 x 10"^'^ erg cm~^ 
(Paper I). This contrast is unprecedented for an AXP. 
Further RXTE observations yielded SGR-like episodes 
of bursts (IWoods et alJl200,^. similar to those seen from 
IE 2259-H586 l|KasDi et al.ll2003j) . A search for an opti- 
cal/IR counterpart detected a fading IR source within 
the Chandra error circle, similar to ones associated with 
other AXPs, confirm ing its identification with XTE J1810 
(|R.ea et a1.ll2004albD . 

In this paper we present the latest results from TAXP 
XTE J1810-197 including new XMM-Newton observa- 
tions extending to 2006 March. These observations allow 
us to characterize the spectral evolution of a TAXP in 
outburst as it returns to quiescence. We show that the 
spectrum has now begun a marked transition back to 
the nominal quiescent state, with a distinct temperature 
and flux evolution. This long-term evolution provides 
the unique possibility of decomposing its fading spectral 
comp onents in a manner unav ailable for other magne- 
tars l|HalDern fc Gotthel j l2005t hereafter Paper II). In 
the following, we used a revised distant to the NS, o?3.3, 
quoted in units of 3.3 kpc and based on the radio pulsar 
DM measurement of .Camilo et al.i (|2006|') . 



2 Long-term spectral and temporal evolution 

Herein we present a total of seven XMM-Newton obser- 
vations of XTE J1810— 197of whic h four have been pre- 
viousl y described in Papers I, II, and lGotthelf fc HalperiJ 
l)2005l hereafter Paper III). The three new data sets were 
obtained using similar observing modes as for previous 
observations and reduced and analyzed in an identical 
manner. A log of these observations is recorded in Ta- 
ble 1 . A full report on the reduction and analysis of these 
data sets will be presented in Gotthelf et al. (in prep.). 

Figure Q presents an up-to-date light curve of TAXP 
XTE J1810-197 derived by adding the XMM-Newton 
flux m easurements to the RXTE results of llbrahim et al.l 
l)2004() . For comparison, the XMM-Newton fluxes were 
extracted assuming a simple power-law spectral model 
fitted in the 2 — 10 keV energy band. The combined data 
points are well fitted by an exponential decay model with 
overall time-constant of 233.5 d, somewhat shorter but 
consistent with the initial RXTE trend, given the ex- 
tended monitoring interval. To allow for a systematic 
offset in the RXTE flux measurements^ , these data were 
rescaled by a factor of 1.42 to match the XMM-Newton 
results. The implied flux range for the initial outburst is 

^ llbrahim et ai.l 120041) assumed a nominal 2 — 10 keV 
conversion factor of 2.27 counts per CPU = 2.4 x 
10"" erg s"^ cm"^. 
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Fig. 1 Long-term light-curve of XTE J1810— 197 using mea- 
surements obtained with RXTE (squares) and XMM-Newton 
(diamond). The RXTE data are from Ibrah im et a l. (2004), 
renormalized to match the XMM-Newton results. The XMM- 
Newton data points are for the observations of Table 1, but 
fitted with a power-law model in the 2 — 10 keV energy band, 
for direct comparison with the RXTE results. The solid line 
is a best combined fit to an exponential decay model (see text 
for parameters). 

F(2 - 10 keV) = (11 - 8) x 10"" erg s^^ cni^^ over 
the interval during which the target was inaccessible to 
RXTE. Clearly the flux in this band has all but returned 
to its pre-outburst value. 

Spectra of AXPs are nominally modeled assuming a 
two-component power-law plus blackbody model. The 
spectral and temporal trends for XTE J1810— 197 al- 
low us to strongly reject this model (see Paper II for a 
comprehensive discussion). Instead, we find that a two- 
temperature blackbody model gives equally acceptable 
spectral fits but is better motivated physically (see §5). 
A summary of spectral results is presented in Table 1 
and the spectra, fitted with the double blackbody model, 
are displayed in Figure |21 (the 2003 Oct observation is 
excluded here for clarity). Although the hot blackbody 
component (BB2; Fig. |2Jl initially dominated the emis- 
sion at low energies, it fades more rapidly than the cooler 
"warm" emission component (BBl; Fig. Evidently 
the spectral fits after 2004 March deviate significantly 
at the lower end of the XMM-Newton energy band, be- 
low 0.7 keV. It is not yet clear if this is an instrumental 
artifact, however, taken at face value there seems to be 
an absorption feature at ~ 350 eV with a Gaussian width 
of cr « 160 eV. The mean equivalent width of this feature 
is 600 eV, but its strength is somewhat time dependent. 
Perhaps as the hotter component contributes less, this 
feature becomes more evident, suggesting an association 
exclusively with the cooler component. Further research 
is needed to quantify this feature and consider its reality, 
or physical implications. 

The spectrum of XTE J1810-197 can be thought of 
as the combined flux from two concentric emitters, "hot 
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Table 1 XMM-Newton Spectral Results for XTE J1810-197 



Parameter 


2UU6 bep 


2OO0 Uct Iz 


2U04 Mar 11 


2004 bep lo 


2U0o Mar lo 


2U0o bep 2U 


2U0D Mar lo 


Expo (ks)" 
iVu (cm-^)'' 


11.5/8.1 


6.9/6.2 


17.0/15.8 


26.5/24.4 


39.8/37.2 


39.5/37.8 


41.7/38.8 


0.65 ± 0.04 


0.65 ± 0.04 


0.65 ± 0.04 


0.65 (fixed) 


0.65 (fixed) 


0.65 (fixed) 


0.65 (fixed) 


kTi (keV) 


0.25 ± 0.02 


0.29 ± 0.04 


0.27 ±0.02 


0.25 ±0.01 


0.22 ±0.01 


0.20 ±0.01 


0.19 ±0.01 


fcTa (keV) 
Ai (cm^) 


0.68 ±0.02 


0.71 ±0.03 


0.70 ±0.01 


0.67 ±0.01 


0.60 ±0.01 


0.52 ±0.01 


0.46 ± 0.02 


5.6 X 10^^ 


2.9 X 10^^ 


3.3 X 10^^ 


4.0 X 10^^ 


4.9 X 10^^ 


6.6 X 10" 


7.2 X 10" 


A2 (cm^) 


2.8 X 10" 


2.2 X 10" 


1.3 X 10" 


8.7 X 10"' 


6.0 X 10"' 


3.7 X 10"' 


3.6 X 10^" 


BBl Flux= 


4.2 X 10"^^ 


5.4 X 10"^^ 


3.5 X 10"" 


2.6 X 10~" 


1.6 X 10"" 


1.0 X 10"" 


7.5 X 10"^^ 


BB2 Flux'= 


3.5 X 10"" 


3.0 X 10~" 


1.8 X 10"" 


1.0 X 10"" 


4.0 X 10"" 


1.3 X 10~" 


6.8 X 10"" 


Total Flux'^ 


3.93 X 10"" 


3.84 X 10"" 


2.13 X 10"" 


1.29 X 10"" 


5.67 X 10"" 


2.35 X 10"" 


1.44 X 10"" 


LBBl(bol)'' 
LBB2(bol)'' 


2.4 X lO^'' 


2.3 X lO^'' 


1.7 X lO^'' 


1.6 X 10^" 


1.2 X lO^'' 


1.0 X 10^^ 


8.6 X 10^^ 


6.3 X lO^'* 


5.7 X 10^" 


3.1 X 10^" 


1.8 X 10^** 


7.9 X 10^^ 


2.8 X 10^^ 


1.7 X 10^^ 




1.1(187) 


1.1(84) 


1.1(194) 


1.2(188) 


1.6(152) 


1.6(80) 


1.6(117) 



Note - Uncertainties in spectral parameters are 90% confidence for two interesting parameters. 
'^EPIC-pn exposure/livetime in units of ks. 



'interstellar hydrogen absorbing column density in units of cm~^. 
■^Absorbed 0.5-10 keV flux in units of erg cm"^ s~^. 

■^Unabsorbed bolometric luminosity in units of erg s~^ assuming a distance of d = 3.3 kpc. 



spots" , whose temperature and size are evolving at dif- 
ferent rates, effectively changing the overall shape of the 
spectrum with time (Paper II). With a set of spectral 
measurements spanning two and a half years, a clear 
trend has emerged. The bolometric luminosities of the 
two components are shown in Figure 13 This reveals a 
cooler component decreasing exponentially in time with 
Ti = 870 d, while the hotter temperature flux declines 
with T2 — 280 d. As might be expected, the shorter 
time constant is very close to the one described by the 
long-term flux above 2 keV (Figure , where the hot- 
ter blackbody component dominates. However, the lat- 
est two data points show that this component has now 
fallen faster than the original exponential. Although the 
bolometric luminosities of the two components are quite 
different at the latest epoch, their measured fluxes are 
now nearly equal in the 0.5 — 8 keV range. We also note 
that, while it is possible to fit an alternative power-law 
decay model to the hot blackbody flux, such a model 
would require a decay index that steepens with time. 

Based on the decay rates and approximate outburst 
time we estimate an initial bolometric luminosity for the 
two spectral components of ii « 3 x lO'^'^dg g erg s~^ 
and L2 ~ 2 X 10^^ 3 erg s~^, respectively. The total 
luminosity at the peak of the outburst is comparable to 
that of a persistent AXP, suggesting that this TAXP is 
anomalously faint. We can now compute the implied flu- 
ences of /i « 2 x 10**^ d| 3 erg and /2 « 4 x 10'*^ (i§ 3 erg, 
for the two components, respectively. The energy of the 
outburst event is many orders-of-magnitude lower than 
that for a typical SGR flare. 

The temperatures and inferred blackbody emitting 
areas have also evolved measurably with time. Figure^ 
displays the time histories of these parameters for each 
spectral component, obtained from the model fits. Based 




2004 2005 2006 2007 

TIME (YEAR) 



Fig. 3 A semi-log plot of the bolometric luminosity of 
XTE J1810— 197 as a function of time (crosses) for each tem- 
perature component of the double blackbody model given in 
Table 1. The light-curves are fit assuming an exponential de- 
cays model (solid lines); the corresponding e- folding times is 
Ti = 280 d and T2 = 870 d for the warm and hot temperature 
components, respectively. The recent data for the hot com- 
ponent now deviates from this model. The fitted quantities 
have been extrapolated to the (la) quiescent range measured 
in Paper I (cross-hatched area). 



on the three latest data points, the temperatures now 
show a definite cooling trend for each spectral compo- 
nent, as suggested by the broken line in Figure^l Prior to 
mid 2004, the temperatures likely remained nearly con- 
stant, after which they both fell at a rate of ~ 22% per 
year (AkTi = -0.051 keV yr'^; AkT2 = -0.15 keV yr'^). 
The size of the effective hot spots (derived from the 
blackbody emission areas) also followed distinct evolu- 
tions. The hot component has been shrinking exponen- 
tially since the initial XMM-Newton observation, while 



4 



E. V. Gotthelf, J. P. Halpern 



2003 SEP 08 





ENERGY (KEV) 




2 5 
ENERGY (KEV) 



(c) 




2004 SEP 18 : 






BB2\^ - 










/ 2 


5 



ENERGY (KEV) 



3 



3 



2005 MAR U 




2 5 
ENERGY (KEV) 



: (e) 


2005 SEP 20 : 




BB1 ) BB2 ^ 







2 5 
ENERGY (KEV) 



(J) 



2006 MAR 12 




ENERGY (KEV) 



Fig. 2 XMM-Newton EPIC pn spectra of XTE J1810~197 from the earliest to tiie latest epochs, in six months intervals 
(the 2003 Oct spectrum is excluded, for clarity). These spectra are shown with the best-fit two-temperature blackbody model 
specified in Table 1. Although the temperatures of the blackbody components have not changed greatly between epochs, the 
flux of the hot component (BB2) has decayed rapidly relative to that of the cooler one (BBl). Also shown are the residuals 
from the best-fit models. The nature of the deviations to the model below 0.7 keV has yet to be determined (see text). 



the warm component has steadily increased in size. A no- 
table exception is the initial data point which deviates 
from the trend for both the areas and temperatures. This 
may be associated with a glitch or rotation instability, as 
suggested by the pulse timing results around this epoch 
(§3). 

The area of the warm component [Ai] Table 1) may 
be reaching a maximum, corresponding to the whole NS 
surface. This would then provide a lower-limit on the 



NS radius of ^ 7.5 da. 3 km, ignoring relativistic effects 
(redshift and light bending). In contrast, the area of the 
hotter component continues to shrink and its contribu- 
tion to the flux is severely diminished. 

As shown in Figure |21 the X-ray luminosity of TAXP 
XTE J1810— 197 has nearly returned to its historic qui- 
escent level. This is also true of the temperature and 
inferred blackbody area. The pre-burst ROSAT spec- 
trum of 1992 March 7 is reasonably well fitted with 
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2004 2004.5 2005 2005.5 2006 
YEAR 

Fig. 4 The time history of the best-fit XMM-Newton 
component temperatures and surface emitting areas for 
XTE J1810-197 based on the double blackbody spectral 
model (Table 1). The evolving areas (Top) indicates that the 
warm component is expanding to cover the NS surface, while 
the hot component is shrinking rapidly. The corresponding 
temperatures {Bottom) are declining at a steady rate of 22% 
per year since mid 2004. Prior to that time the temperature 
of the hotter component was essentially constant. This trend 
for the warm component is somewhat ambiguous. 



a single blackbody of kT = 0.18 ± 0.02 keV covering 
5.2 X 10^2 cin2 and LBB(bol) = 5.6x10^3^2 3 ergs"! 
(Paper I). The latest XMM-Newton spectrum that over- 
laps the ROSAT band is again dominated by the cooler 
blackbody component. Furthermore, the current pulsed 
fraction is consistent with the quiescent upper-limits. 
Thus the double blackbody model is sufhcient to de- 
scribe the observed properties, and may not require a 
third "quiescent" component to explain the pre-outburst 
observations. XMM-Newton data from the next observ- 
ing window (2006 September) is likely to be consistent 
with the ROSAT result in the soft energy band, effec- 
tively defining a return to the quiescence state. 





















































XTE Jiai0-197 
SPIN-DOWN HISTORY 




dP/dt = (1.26±0.04) X 



2004 2004.5 2005 2005.5 2006 
YEAR 

Fig. 5 The spin-down history of XTE J1810-197 pulse pe- 
riod as measured using the XMM-Newton data of Table 1. 
The initially erratic spin-down appears to be settling down - 
the last 5 data points are well-modeled by a linear spin-down 
model {solid line); in contrast, a poor fit is obtained using all 
points {dotted line). 



3 Spin-down Evolution of XTE J1810-197 

The barycentered pulse period of XTE J1810— 197 mea- 
sured at each XMM-Newton observing epoch is shown 
in Figure \^ These data points were derived following 
the method outlined in Paper I. Prior to epoch 2004, 
the spin-down rate derived from the RXTE and initial 
XMM-Newton measurements were highl y erratic, rang- 
ing fr om P = (0.8 - 2.2) x 10"" s s'^ l|Ibrahim et alJ 
120041: Paper III). This temporal behavior is likely asso- 
ciated with the outburst event. The latest period mea- 
surements suggest that the spin-down rate has settled 
down somewhat. A linear model fit to the last 5 pe- 
riod measurements yields a period derivative of P = 
1.26 ± 0.04 X 10"" s s-i with e. xl = 1-2 for 3 DoF, 
corresponding to a null hypothesis probability of 0.27. 
Compared to the radi o P = 1.016 ± 0.001 x 10~" s s"^ 
of 17 March - 7 May ijCamilo et al J 12006). there is evi- 
dently still timing noise in the pulsar's spin-down. The 
revised X-ray rate implies a nominal characteristic age 
Tc = 7.0 kyr, surface magnetic field Bs = 2.7 x lO^'' G, 
and spin-down power E = 2.9 x lO^^ erg s^^, values 
typical for a magnetar. 

The energy-dependent modulation of the pulsar has 
evolved noticeably over time with the declining flux. Fig- 
ure El displays the pulse profiles at each XMM-Newton 
epoch, folded at the best determined periods. The data 
has been broken up into six energy bands, with the cor- 
responding background in each band subtracted. The 
pulsed fractions, defined as the pulsed emission divided 
by the total flux, are indicated on the corner of each 
panel of the figure. These profiles are aligned so that 
phase zero is the same at each epoch, for ease of com- 
parison. In any case, within an observation, the phase 
zero does not change with energy. This suggests a geo- 
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Instead we conclude that either the two spectral compo- 
nents contain an admixture of the two shapes or there 
is a third, unmodeled spectral component present. This 
is a direction of active research, to consider the correct 
characterization of the phase dependent spectral contri- 
bution. 



4 Radio observations of XTE J1810-197 



metric interpretation for the modulation, such as from 
localized emission on the rotating NS star surface. The 
pulsed fraction has generally decreased with time, most 
notably at X-ray energies below E < 2 keV; at higher 
energies the trend is less clear due to the large uncer- 
tainties derived for those pulsed fractions. Furthermore, 
at all epochs, the pulsed fraction clearly increases with 
energy, an equally interesting result discussed below. 

The evolving spectral shape and pulse profiles of TAXP 
XTE J1810— 197 provide an important and unique (so 
far) diagnostic of the emission geometry, and ultimately, 
the emission mechanism(s) of NSs. Because the spectrum 
is well modeled by two blackbody components, and the 
phase alignment of the pulse profiles are energy inde- 
pendent, it is most natural to consider emission from 
two concentric regions on the NS surface. The hot com- 
ponent is associated with a smaller hot-spot, while the 
blackbody model predicts a larger annulus for the warm 
component. The complete collection of observed pulsed 
fractions is consistent with this model. The smaller, hot- 
ter spot always dominates the spectrum above 2 keV and 
offers a natural explanation for the higher modulation at 
these energies. However, at lower energies, this contribu- 
tion gradually fades relative to cooler emission (which is 
less modulated), contributing less and less to the pulsed 
emission over time. 

We note in passing that a fitted power-law spectral 
component would have to dominate the soft (< 2 keV) X- 
rays at all epochs, while varying in its contribution to the 
hard (> 2 keV) X-rays. Such a power law would drive an 
evolution of the pulse shapes that is opposite of what is 
observed. Thus, we find that the detailed evolution of the 
X-ray emission from XTE J1810— 197 further supports 
the assumption of a purely thermal spectral model, and 
leaves no evidence of a steep power law as is commonly 
fitted to individual observations of AXPs. 

Since the first observation of XTE J1810-197 it was 
apparent that the broad-band pulse shape is not a simple 
sine function; the pulse peak is relatively sharp with a 
broader inter-pulse trough (Paper I). This effect is more 
pronounced at higher energy, where the profile is nearly 
triangular in shape. This suggests that the pulse profiles 
can be decomposed into a triangle function and a sinu- 
soidal function, a model that was explored for the earlier 
data sets and detailed in Paper II & III. This model con- 
tinues to be appropriate for the new data and is used to 
extract unbiased pulsed fraction measurements for the 
profiles shown in Figure 

Given the success in modeling the pulse profiles with 
these functions it is natural to associate the two pulsed 
components uniquely with the two spectral components, 
i.e., the triangle shape for the hotter blackbody compo- 
nent and the sinusoidal shape for the cooler one. How- 
ever, we find that we can not model all the profiles in 
a consistent manner with just a simple superposition of 
these temporal components, based on the implied fiu x 
ratio in each energy band ijCrotthelf fc HalnernlEinl . 



Considering that the absence of radio emission is a defin- 
ing characteristic of AXPs, the serendipitous radio detec- 
tion from XTE J1810— 197 came as a great surprise. A 
chance search of VLA data taken about a year after out- 
burst reveals an unresolved point source with a flux of 
4.5 mJy at 1.43 GHz, located at the precise coordinates 
of the TAXP iHaloern et al.l200fibi) . Other archival VLA 
data obtained at various frequencies provide upper-limits 
before and after outburst. A follow-up VLA observation 
in 2006 March yielded a flux of 12.9 mJy at 1.43 GHz. 
Together, these results indicate highly erratic radio emis- 
sion. A second surprise came with the detection around 
that time o f pulsed radio emission at the X-ray period 
(ICamilo et a l. 200gJ. This search, at = 1.4 GHz using 
the Parkes radio observatory, revealed a narro w, bright 
pulse with a high degree of linearly polarization ijCamilo et alJ 
l2006j) . A series of multi- frequency measurements in 2006 
April-May confirms the erratic flux behavior and pro- 
vides an unusual spectral slope of a ^ — 0.5, where 
Si, oc i^" (cf. a ^ —1.6 for a typical radio pulsar). A re- 
vised distance to the pulsar of 3.3 kpc is inferred from the 
pulsar's dispersion. A search for pulsations in archival 
data acquired in 1997/1998 produced a null result and 
argues against significant emission prior to the outburst 
event. A summary of these radio observations is shown 
in Figured 



5 Emission Geometry: Models & Theory 

Although transient AXPs such as XTE J1810-197 are 
relatively rare, their short active duty cycle suggests the 
existence of a larger population of unrecognized young 
NSs. XTE J1810— 197 provides a unique window into 
this population, with prior measurements in the quies- 
cent state and detailed observations during its active, 
pulsed phase. Ultimately, we hope to use this TAXP to 
give insight into the emission mechanisms of magnetars, 
in general. With the discovery of pulsed radio fiux, this 
may carry over to interpreting emission mechanism of ra- 
dio pulsars, as well. Below we now summarize our initial 
attempts to interpret the XTE J1810— 197 results in the 
context of a physical model. First we outline the basic 
arguments for the two-temperature blackbody spectral 
model as a plausible alternative to the nominal power- 
law plus blackbody model. 
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Fig. 6 Energy-dependent pulse profiles of XTE J1810— 197 obtained with the XMM-Newton EPIC pn detector for all seven epochs of Table 1. The background 
for each profile has been subtracted and phase zero is aligned for each epoch, for comparison. The phase of the peak is seen to be energy independent. The pulsed 
fraction at low X-ray energies has decreased with time, while remaining essentially unchanged at high energy. Also shown is the best fit to the two-component 
model for the pulse profile {solid line) described in the text (see §2.1). 
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Fig. 7 A set of archival and dedicated radio observations of 
XTE J1810— 197 measured at three frequencies. Prior to out- 
burst (denoted by vertical lines) no radio emission is found 
corresponding to the plotted upper-limits. Following the X- 
ray outburst, large variations in flux are inferred from sev- 
eral measurements and upper-limits. Pulsed radio emission 
is also recorded at the X-ray period in the post-outburst 
era following the 2006 discovery. A search of archival data 
prior to the outburst places an upper-limit of ~ 0.4 mjy at 
1.4 GHz. For clarity, multi-frequency measurements coincid- 
ing with the 1.4 GHz puls ed emission data points are not 
shown (from lHalpern et anr2005bi: iCamilo et al.ii200Q') . 



includes general relativistic effects (redshift and light 
bending) . Our goal is to match spectrum and energy de- 
pendent pulse shape in order to determine viewing ge- 
ometry, distance, and NS size. Preliminary work is able 
to reproduce the pulsed fraction to a reasonable degree 
but the pulse shape remains elusive. 

A framework for a theoretical interpretation of the 
emission from XTE J 1810— 19 7 is suggested by the mag- 
netar coronal model of lBeloborodov fc Thompson. (.2006) . 
According to this model, the large outburst was gener- 
ated by a starquake that resulted in a transition to an 
active coronal state that caused energy to be stored in 
the twisted B-field of the coronal loop. Particles (mostly 
e+e~) are accelerated within this loop and impact the 
NS surface with GeV energy. This heats up the loop 
footprint resulting in the observed hot spot emission. 
The decay timescale of the coronal loop, and thus the 
hotter component of the double blackbody luminosity, 
is of order of a few years. The decay rate is determined 
by ohmic dissipation of current in the excited loop. A 
cooler component likely arises from deep crustal heat- 
ing associated with the initial outburst and possibly ear- 
lier ones. The features of this model are in general ac- 
cordance with the observational properties and inferred 
model for XTE J1810-197. 



While both models give reasonable fits to the spec- 
trum, it's important to note that the power-law compo- 
nent is used to model the low energy residuals and not 
any high- energy tail. Without a low energy cut-off of this 
component, it is not possible to connect it with the ob- 
served IR emission without an energy catastrophe. As de- 
tailed in Paper II, invoking synchrotron self-absorption is 
excluded by radius/magnetic-field inconsistently. In con- 
trast, the extrapolated spectrum of the warm blackbody 
component does not exceed the measured IR flux. Fur- 
thermore, currently there is no acceptable physical model 
to anchor a power-law component. As discussed earlier, 
the two components of the double blackbody model can 
be associated with a pair of hot spots on the surface of 
the star. The warm temperature component, covering a 
large fraction of the NS surface, is consistent with the 
decreased pulsed fraction at lower energies. The hotter 
blackbody component is consistent with a smaller emis- 
sion area and greater modulation at higher energies. A 
comprehensive discussion of this issue can be found in 
paper II. 

Determining the emission geometry on the NS is of 
great interest. The natural conclusion of applying the 
double blackbody model to the time resolved spectra 
and pulse profiles is that concentric hot regions give rise 
to the observed modulation. We consider a model of the 
phase-resolved spectrum taking into account the viewing 
geometry and offset of the hot spots from the rotation 
axis (Perna & Gotthelf 2006, in prep.). This is bas ed on 
the NS emission model given in lPerna et al.l l)200lj) that 
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